home *** CD-ROM | disk | FTP | other *** search
-
-
-
- ____BBBBAAAALLLLAAAANNNNCCCC((((3333FFFF)))) ____BBBBAAAALLLLAAAANNNNCCCC((((3333FFFF))))
-
-
-
- NNNNAAAAMMMMEEEE
- BALANC, SBALANC - EISPACK routine. This subroutine balances a REAL
- matrix and isolates eigenvalues whenever possible.
-
-
- SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
- ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee bbbbaaaallllaaaannnncccc((((nnnnmmmm,,,, nnnn,,,, aaaa,,,, lllloooowwww,,,, iiiigggghhhh,,,, ssssccccaaaalllleeee))))
- iiiinnnntttteeeeggggeeeerrrr nnnnmmmm,,,, nnnn,,,, aaaa,,,, lllloooowwww,,,, iiiigggghhhh
- ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn aaaa((((nnnnmmmm,,,,nnnn)))),,,, ssssccccaaaalllleeee((((nnnn))))
-
- ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssbbbbaaaallllaaaannnncccc((((nnnnmmmm,,,, nnnn,,,, aaaa,,,, lllloooowwww,,,, iiiigggghhhh,,,, ssssccccaaaalllleeee))))
- iiiinnnntttteeeeggggeeeerrrr nnnnmmmm,,,, nnnn,,,, aaaa,,,, lllloooowwww,,,, iiiigggghhhh
- rrrreeeeaaaallll aaaa((((nnnnmmmm,,,,nnnn)))),,,, ssssccccaaaalllleeee((((nnnn))))
-
-
-
- DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
- On INPUT
-
- NNNNMMMM must be set to the row dimension of two-dimensional array parameters
- as declared in the calling program dimension statement.
-
- NNNN is the order of the matrix.
-
- AAAA contains the input matrix to be balanced. On OUTPUT
-
- AAAA contains the balanced matrix.
-
- LLLLOOOOWWWW and IGH are two integers such that A(I,J) is equal to zero if
- (1) I is greater than J and
- (2) J=1,...,LOW-1 or I=IGH+1,...,N.
-
- SSSSCCCCAAAALLLLEEEE contains information determining the
- permutations and scaling factors used. Suppose that the principal
- submatrix in rows LOW through IGH has been balanced, that P(J) denotes
- the index interchanged with J during the permutation step, and that the
- elements of the diagonal matrix used are denoted by D(I,J). Then
-
- SSSSCCCCAAAALLLLEEEE (J) = P(J), for J = 1,...,LOW-1
- = D(J,J), J = LOW,...,IGH
- = P(J) J = IGH+1,...,N. The order in which the interchanges
- are made is N to IGH+1, then 1 TO LOW-1. Note that 1 is returned for IGH
- if IGH is zero formally. The ALGOL procedure EXC contained in BALANCE
- appears in
-
- BBBBAAAALLLLAAAANNNNCCCC in line. (Note that the ALGOL roles of identifiers K,L have been
- reversed.) Questions and comments should be directed to B. S. Garbow,
- Applied Mathematics Division, ARGONNE NATIONAL LABORATORY
-
-
-
-
-
-
-
- PPPPaaaaggggeeee 1111
-
-
-
-